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Abstract: 

A previous calculation of electroweak 0(a) corrections to deep- inelastic neutrino scat- 
tering, as e.g. measured by NuTeV and NOMAD, is supplemented by higher-order effects. 
In detail, we take into account universal two-loop effects from Aa and Ap as well as higher- 
order final-state photon radiation off muons in the structure function approach. Moreover, 
we make use of the recently released (9(a)-improved parton distributions MRST2004QED 
and identify the relevant QED factorization scheme, which is DIS like. As a technical 
byproduct, we describe slicing and subtraction techniques for an efficient calculation of 
a new type of real corrections that are induced by the generated photon distribution. A 
numerical discussion of the higher-order effects suggests that the remaining theoretical un- 
certainty from unknown electroweak corrections is dominated by non-universal two-loop 
effects and is of the order 0.0003 when translated into a shift in sin 2 #w — 1 — M^/M^. 
The 0(a) corrections implicitly included in the parton distributions lead to a shift of 
about 0.0004. 
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1 Introduction 



In 2002 the NuTeV collaboration has deduced the on-shell weak mixing angle, 
sin 2 6* w = 1 — M^/M|, from the ratio of neutral- to charged-current cross-sections mea- 
sured in deep-inelastic neutrino scattering [1]. The NuTeV result on sin 2 # w is about 3a 
at variance with the overall fit [2] of the Standard Model (SM) to precision observables, 
a fact that is known as "the NuTeV anomaly" . In spite of many attempts to understand 
or to resolve this anomaly in or beyond [3-7] the SM, the issue is still not settled. Inside 
the SM, uncertainties originating from parton distributions [4,5,7-10] (isospin violation, 
strange-sea asymmetry), from nuclear effects [4,7,9,11], and from QCD [3,9,12] and elec- 
troweak [13-16] radiative corrections significantly influence the NuTeV result. In view of 
the recent progress on these subjects, a careful reanalysis of the NuTeV result, including 
its error assessment, would be welcome. In this context it is very interesting that the 
NOMAD collaboration plans to extract sin 2 # w from their data on neutrino scattering 
[17]- 

In our previous work [15], we performed a calculation of the electroweak corrections of 
0(a) to neutral- and charged-current deep-inelastic neutrino scattering and compared our 
results with the older calculation [14] which was used in the NuTeV data analysis. Unfor- 
tunately, we could not fully reproduce the results of Ref. [14], which was in part certainly 
due to the incompletely known setup and input. However, a comparison with a more 
recent recalculation [16], which is in line with Ref. [14], suggests that the still remaining 
differences are due to different kinematical approximations made in Refs. [14,16] compared 
to our calculation. While our calculation [15] is a pure 0(a) calculation without further 
approximations, in Refs. [14,16] the momentum transfer of the W and Z bosons was ne- 
glected, i.e. the approximation of a four-fermion contact interaction was made. Moreover, 
we concluded in Ref. [15] that the estimate of remaining theoretical uncertainties from 
unknown electroweak higher orders was too optimistic. 

In this paper, we refine our previous calculation [15] of the electroweak corrections to 
deep-inelastic neutrino scattering by including leading effects beyond 0(a), in order to 
further reduce the theoretical uncertainty in the electroweak sector. In detail, after a short 
summary of our conventions in Section 2 we work out the leading universal corrections 
originating from the running of the electromagnetic coupling and from the p parameter in 
various input-parameter schemes in Section 3, in order to reduce the scheme dependence. 
Moreover, in Section 4 we include the leading logarithmic effects beyond 0(a) that result 
from collinear final-state radiation (FSR) in the structure-function approach. Finally, we 
make use of the recently released (9(a)-improved parton distribution functions (PDFs) 
of the MRST collaboration [10], which requires two kind of refinements, as explained in 
Section 5.1. Firstly, the factorization scheme for the 0(a) corrections has to be chosen 
in a way that is consistent with the treatment of 0(a) corrections used in the fit of 
the PDFs. Secondly, the 0(a) driven evolution of the PDFs leads to an induced photon 
distribution, so that a new type of real corrections with a photon in the initial state has to 
be considered. For a proper technical treatment of the photon-induced processes we have 
worked out different variants of phase-space slicing and dipole subtraction, as described 
in the Appendix in detail. Section 6 contains our numerical results, which comprise a 
discussion of the various higher-order corrections both to integrated and differential cross 
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sections. Moreover, we derive an estimate of the remaining theoretical uncertainties from 
still missing electroweak effects there. Our conclusions are given in Section 7. 



2 Kinematics and conventions 

In this section we briefly repeat our conventions of Ref. [15] as far as relevant in the 
sequel. We consider the neutral-current (NC) and charged-current (CC) parton processes 

NC: i/p(pi) + q(p q ) -> Vpih) + q(k q ), q = u, d, s, c, u, d, s, c, (2.1) 
CC: Vp(pi) + q(p q ) -> n'ih) + q\k q ), q = d, s, u, c, g' = u,c,d,s, (2.2) 

and the processes with all particles replaced by their antiparticles. The assignment of 
momenta is indicated in parentheses. At lowest order the former processes proceed via 
Z-boson exchange, the latter via W-boson exchange. CP symmetry implies that the 
(parton) cross sections for z/^-quark and z/^-antiquark scattering are equal to the ones of 
z/^-antiquark and z/^-quark scattering, respectively. 

In the discussion of numerical results, we concentrate on the NC to CC cross-section 
ratio [18] 

RV = a^N^u,X) 

aZd^N -> p-X)' V ' 

but our calculation is applicable to any combination of NC and CC cross section involving 
incoming neutrinos or anti-neutrinos, such as the ratio 

" -> ijl-X) - oZcfrN - iJL+XY { ■ ] 

as proposed by Paschos and Wolfenstein [19]. 

Throughout the calculation, we neglect the muon mass and the parton masses when- 
ever possible. The quark masses are only kept as regulators for mass singularities, which 
appear as mass logarithms in the (photonic) corrections; the quark-mass logarithms are 
absorbed into the PDFs as described in Section 5. The muon mass is not only required 
as regulator for collinear final-state radiation but also for a proper description of forward 
scattering in the (loop-induced) 7^z/ M vertex, as described in Ref. [15] in detail. 

The momentum of the incoming (anti-)quark q is related to the total nucleon momen- 
tum P N by the usual scaling relation 

< = xP&, (2.5) 

which holds in the centre-of-mass (CM) frame of the partonic system where the nucleon 
mass Mtv is negligible. The variable x is the usual momentum fraction, restricted by 
< x < 1. Since deep-inelastic neutrino scattering usually is accessed in fixed-target 
experiments (as NuTeV), we take the energy Ejj AB of the incoming (anti-)neutrino beam 
to define the incoming momentum pi in the rest frame of the nucleon, called LAB in the 
following. Thus, the squared partonic CM energy is given by 

s=( P i+ Pq ) 2 = 2xM N E^ AB , (2.6) 
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up to terms of higher order in the nucleon mass. Moreover, we define the Bjorken variable 



-(Pi - h? 
2P N ( Pl - ki) 



-(Pi 



k,) 2 



had+phot 



and the energy ratio y via 



y 



pLAB 

-^had+phot 

^LAB • 



(2.7) 



(2.8) 



where -Ehad+phot i s the energy (in the LAB frame) deposited by hadrons and photons in 
the detector. In lowest order xbj and x coincide, and -Ehad+phot is identical with the energy 
of the outgoing quark; in the presence of photon (or gluon) radiation this is not the case 
anymore. 



3 Universal electroweak corrections beyond O(ct) 

The integrated CC and NC cross sections at parton level, including 0(a) corrections, 
can be written as 



°"nc,i 



IPS 



IPS 
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= a 



NC,0 
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(3.1) 



where o"cc,o and a^c o denote the partonic lowest-order cross sections, and r = ± refers to 
the chirality of the quark line. The input-parameter scheme (IPS) is indicated by an upper 
index which takes the value a(0), a(M z ), or for the respective scheme. The precise 
definition of these IPS is given in Section 4.2 of Ref. [15]. Roughly speaking, the whole 
SM input consists of all particle masses together with the chosen a for the electromag- 
netic coupling. The a(0)-scheme corresponds to the "complete on-shell renormalization 
scheme" as described, e.g., in Refs. [20,21]. The relative 0(a) corrections <5cc,i and S^ c 1 
are proportional to the value of a in the corresponding IPS, where q;g m = V^G^M^s^/ir 
in the G M -scheme. Apart from this factor the relative corrections differ by constant terms, 
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Here 
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(3.2) 



(3.3) 



describes the running of the electromagnetic coupling from the scale Q = to Q = M z , 
induced by the light fermions / (all but the top quark) with colour factor Np electric 
charge Qf, and mass rrif. The quantity IT n {Q 2 ) is the photonic vacuum polarization in 
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the convention of Ref. [20]. The quantity Ar 1 is the pure one-loop correction to muon 
decay [22], as defined by 
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1 - 



Ml, 



7ra(0) 

V2G, 



[1 + An(a(0), M w , M z , M H , m/)] 



(3.4) 



Note that Ari implicitly contains large contributions from Aa(Mz) 
loop) correction (c 2 v /s 2 v )Ap 1 ~ 3% induced by the p-parameter, 



6% and the (one- 



An = Aa(M z ) - -fA Pl W +Ar, 



a(0) 



(3.5) 



where 
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(3.6) 



The remainder Ar rem comprises the non-universal one-loop corrections to muon decay and 
is of the size of ~ 1%. 

The leading universal two-loop effects induced by Aa(Mz) and by corrections Ap in the 
parameter p = 1/(1 — Ap) can be included as described in Ref. [23]. The procedure is as 
follows. Starting from lowest-order cross sections in the a(0)-scheme, all orders Aa(Mz)™ 
can be effectively included upon substituting a(0) — > a(M z ) = a(0)/[l — Aa(M z )] as 
electromagnetic coupling. The effects from Ap can be absorbed by replacing and 
by appropriate modifications 



4 + clAp, 



l -si 



(3.7) 



This recipe is correct up to 0(Ap 2 ). As shown in Ref. [23], the introduction of a(M z ) 
and 4> 4 also correctly reproduces the correct terms of C(Ao;(Mz) Ap) in processes with 
four light external fermions. Note that in C(Ap 2 ) both one- and two-loop corrections to 
Ap become relevant; explicitly we use the result 



Ap = 3x t [l + p (2) {Ml/ml) x t 



2«s 
9vr 



y+s) 



3x t 



16tt 2 



G, 



Ap! , (3.8) 



with the function p^ given in Eq. (12) of Ref. [24]. In order to avoid double-counting, 
the one-loop effects of Aa(M z ) and Ap have to be subtracted from the full one-loop 

corrections 5\ ^ ^ introduced above. 

The above recipe can be easily modified to the other two IPS. In the a (M z ) -scheme, 
the replacement a(0) — > a(Mz) has already been done, so that the procedure reduces to 
the substitutions s w — > s w and c w — > c w . In the G^-scheme, one has to recall that «g 



effectively involves a factor a(Mz)s^, so that the recipe reads ac 
and c w — > c w ; as a net effect the combination a Gfi /4 is not modified. 

The lowest-order CC cross section is proportional to a 2 / 's*. Thus, using the above 
recipes we find for the corrected cross sections, which include the full 0(a) corrections as 
well as two-loop improvements from Aa = Aa(Mz) and Ap, the results 
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For the a(0)- and a (M z ) -schemes the two given expressions are equal up to non-leading 
two-loop and leading three-loop terms. Equation (3.11) reflects the well-known fact that 
the G M -scheme, which employs G^, M\y, and Mz as independent parameters, effectively 
absorbs the leading effects from Aa(Mz) and Ap in the coupling of W bosons to fermions. 

The lowest-order NC cross section is proportional to tt 2 (fi , g g ^) 2 /(c w s w ) 2 , where the 
coupling of the Z boson to a quark q with chirality r = ± is given by 

r3 

(3.12) 
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9qqZ = Qq + T -, 



with 1^ — ±| denoting the third component of the weak isospin of a left-handed quark. 
In order to absorb corrections from Ap, we define the modified coupling factors 
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The C(a)-corrected NC cross sections, including the leading two-loop improvements, thus, 
read 



a 



NC,l+LI/ 



a(0) 



= a 



NC,0 



a(0) 



( 3jqZ 



C W 5 V 



\-9qqZ C w S w 



1 



[1-Aaf 



+ "NC,1 



«(0) 



2Aa 



2A Pl 



a(0) 



~ 3° r 

9qqZ S w 



; NC,0 



a(0) 



"(0) / jt a(0) , \ / 

1 + <W -2\A Pl -ApJ M - 
+ 3 Aa 2 + 4 Aa Ap f 1 - S T . N 



c J 3 
w g ,5 

r 3 r_ 
9qqZ S w 



9qqZ S v 



+ Ap 2 3 + 



20,(1-34)^ 



T „5 



NC,1+LL 2 



a(M z ) 



T NC,0 



(3.14) 



/ 9 T qqZ C ™ S ™ \ 2 £ T 
_^-9qqZC™S w J 



a(M z ) 



'NC,0 



+ 



-2Api 



a(M z ) 



1 ~~ ~ 3" ^ 



a(M z ) ( 

1 + <W| -2(^Ap 
+ Ap 2 (s 




1 ' ^ e u 



2c w (l - 3s*)/ 3 



(3.15) 



^NC.l+LL 2 



= a 



NC,0 



= cr 



NC,0 



+ 



G u 



G M 2 



1 + ^0,1 -72" A Pl "Ap 1 




2c w (l + 4) J 3 

a 



9qqZ S v 



(3.16) 



4 Higher-order final-state radiation off muons 

In the CC reaction, where a muon appears in the final state, the emission of photons 
collinear to the outgoing muon leads to corrections that are enhanced by large logarithms 
of the form aln(m 2 /Q 2 ) with Q denoting a scale characteristic for the process. The 
KLN theorem [25] guarantees that these logarithms cancel if photons collinear to the 
muon are treated fully inclusively. Since we, however, apply a phase-space cut on the 
energy of the outgoing muon (or equivalently on the remaining total hadronic+photonic 
energy), contributions enhanced by these logarithms survive. The first-order logarithm 
aln(m 2 i /Q 2 ) is, of course, implicitly contained in the full 0(a) correction, so that Q is 
unambiguously fixed in this order. However, it is desirable to control the logarithmically 
enhanced corrections beyond 0{a). This can be done in the so-called structure function 
approach [26], where these logarithms are derived from the universal factorization of the 
related mass singularity. The incorporation of the mass-singular logarithms takes the 
form of a convolution integral over the lowest-order cross section <tcc,o, 



cccllfsr = f da C c,o(pi,P q ] k, kg) J Q 



dzr^(z,g 2 )e cut (^) 



(4-1) 



where the step function O cut is equal to 1 if the event passes the cut on the rescaled muon 
momentum zk\ and otherwise. The variable z is the momentum fraction describing the 
muon energy loss by collinear photon emission. This treatment, thus, corresponds to the 
experimental situation in a fixed-target experiment where photons in the final state are 
not distinguished from the hadrons that shower in the detector, but the muon is analyzed 
as an isolated particle. For the structure function T^(z, Q 2 ) we take into account terms 
up to 0(a 3 ) improved by the well-known exponentiation of the soft-photonic parts 
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with 7s and r(y) denoting Euler's constant and the Gamma function, respectively, and 
the logarithms 
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The parts solely proportional to a power of correspond to collinear (multi-)photon 
emission off the muon, the exponential factor describing resummed soft-photonic effects. 
The non-logarithmic term "—1" in (3^ accounts for a non-singular universal soft-photonic 
correction. The part proportional to /5 M /5 e corresponds to collinear one-photon emission 
followed by a collinear splitting 7 — > e + e~. Treating this contribution as part of T^(z, Q 2 ) 
means that the event selection is based on a muon momentum zki also in this case, i.e. 
that the e + e~ are considered as part of the hadronic and photonic shower in the detector. 
In this context one should also worry about the possibility of fi + [i~ pairs in the final state 
resulting from photon splitting. If one assumes that such muons escape the detector, i.e. 
are not part of the hadronic shower, the corresponding logarithmic contributions cancel, 
since the full z range is integrated over. 

Technically we add the cross section (4.1) to the one-loop result and subtract the 
lowest-order and one-loop contributions da cc LL i FSR already contained within this for- 
mula, 

ct cc,ll 1 fsr = / ^cc,o(pi,Pq',ki,k q ) J dz [5(1 - z) + T^\z,Q 2 )] Q cut (zki), (4.4) 

in order to avoid double counting. The one-loop contribution to the structure function 
reads 



tLL,1 

mm 



z : Q 2 




(4.5) 



More precisely, we adapt the value of a in Y l ^' 1 {z,Q 2 ) to the chosen IPS, so that the 
aln(m 2 /Q 2 ) contribution of the 0(a) correction is subtracted exactly. Thus, the proce- 
dure of adding higher-order FSR effectively changes also the value of a in the a ln(m 2 /<5 2 ) 
term to a(0) which is the appropriate coupling for real-photonic effects. 
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Finally, note that the uncertainty that is connected with the choice of Q 2 enters in 
0(a 2 ), since all 0(a) corrections, including constant terms, are taken into account. As 
default we choose the value 

Q = max^v^, 2m M } (4.6) 

with £ = 1. In order to quantify the scale uncertainty, we vary £ between 0.3 and 3. The 
lower limit 2771^ accounts for the fact that the muon energy \fsj1 in the CM frame cannot 
drop below m M . 



5 Parton distribution functions 



5.1 Inclusion of O(ct) corrections 

The inclusion of 0(a) corrections to hadronic cross sections, which is described in the 
following in detail, conceptually proceeds along the same lines as the incorporation of 
next-to-leading (NLO) QCD corrections. 

The (9(o;)-corrected parton cross sections contain mass singularities of the form 
a\n(m q ), which are due to collinear photon radiation off the initial-state quarks or due 
to a collinear splitting 7 — > qq for initial-state photons. In complete analogy to factor- 
ization in NLO QCD calculations, we absorb these collinear singularities into the quark 
and photon distributions. For processes that involve only quarks and/or antiquarks in 
the initial state in lowest order, the factorization is achieved by replacing the (anti-)quark 
distribution q(x) according to 
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where M is the factorization scale, Cij(z) are the so-called coefficient functions, and the 
splitting functions P%j(z) are defined as 

1 + z 2 
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(5.2) 



The replacement (5.1) defines the same finite parts in the O(a) correction, i.e. the same co- 
efficient functions, as the usual D-dimensional regularization for exactly massless partons, 
where the ln(m g ) terms appear as 1/(D — 4) poles. We have derived this correspondence 
by explicit calculation and found agreement with the results of Ref. [27]. The actual form 
of the coefficient functions defines the finite parts of the 0(a) corrections and, thus, the 
factorization scheme. Following standard definitions of QCD, we distinguish the MS and 
DIS-like schemes which are formally defined by 
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The MS scheme is motivated by formal simplicity, because it merely rearranges the UV- 
divergent terms (plus some trivial constants) as defined in dimensional regularization. 
The DIS-like scheme is defined in such a way that the DIS structure function F 2 does not 
receive any corrections; in other words, the radiative corrections to electron-proton DIS 
are implicitly contained in the PDFs in this case. Whatever scheme has been adopted in 
the extraction of PDFs from experimental data, the same scheme has to be used when 
predictions for other experiments are made using these PDFs. 

The absorption of the collinear singularities of 0(a) into PDFs requires the inclusion 
of the corresponding 0(a) corrections into the Dokshitzer-Gribov-Lipatov-Altarelli-Pa- 
risi (DGLAP) evolution of these distributions and into their fit to experimental data. In 
the past, when the 0(a) effects were not yet included in the PDF fit to data, the inclusion 
of the 0(a) corrections to the DGLAP evolution showed [28] that the impact of these 
corrections on the quark distributions in the MS factorization scheme should be smaller 
than remaining QCD uncertainties. Therefore, 0(a) corrections to hadronic processes 
were usually included in the MS factorization scheme using the existing PDFs, which 
ignored 0(a) effects. This was also done in our previous calculation [15]. 

Recently the MRST collaboration released the first set of PDFs [10], called 
"MRST2004QED" , that consistently include 0(a) corrections. The 0(a) effects induce 
two important modifications in the PDFs. Firstly, the difference in the electric charges 
of the up- and down-type quarks leads to a violation of the isospin symmetry between 
proton and neutron PDFs. It was frequently pointed out in the literature [4,5,7-10] that 
this isospin violation, i.e. v?(x) ^ d n (x) etc., is of particular importance when analyzing 
ratios of neutral-current to charged- current cross sections in neutrino DIS. Secondly, the 
factorization of collinear singularities from photon emission leads to photon distribution 
functions both for the proton and the neutron. The corresponding photon-induced con- 
tributions to hadronic cross sections have not yet been considered in the literature (since 
the photon PDF was not yet available), but have been assumed to be small, because of an 
implicit suppression factor a in the radiatively induced photon PDF. In our discussion of 
numerical results below we illustrate the size of the 0(a) corrections to the PDFs in neu- 
trino DIS cross sections by switching these corrections on and off. Moreover, we calculate 
and discuss the size of the photon-induced contributions to the cross sections explicitly. 
Some technical details of this part of the calculation are given in the Appendix. 

The authors of the MRST2004QED set of PDFs did not explicitly state in Ref. [10] 
which scheme is relevant for the factorization of the 0(a) corrections. They themselves 
did not include any 0(a) corrections, neither in the underlying data set nor in the used 
theoretical predictions, and used the same data as in their standard fit, implicitly assuming 
that no 0(a) corrections corresponding to photon emission off incoming quarks had been 
applied in the data analysis. For the F 2 data analysis at HERA this assumption was 
confirmed to us by members of the HI and ZEUS collaborations and by theorists involved 
in the analysis [29]. Therefore, we conclude that the 0(a) corrections induced by photon 
emission off incoming quarks implicitly went into the PDFs in the fit, which thus precisely 
corresponds to the DIS-like factorization. We use this scheme in the numerical analysis 
below. We note, however, that the MRST2004QED PDFs are defined in NLO QCD in 
the MS factorization scheme. 
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5.2 Nuclear structure functions 

In reasonable approximation, nuclear structure functions can be constructed from 
the PDFs for protons and neutrons upon appropriately weighting the distribution for a 
specific parton according to the probability of finding a proton or neutron in the target. 
The NuTeV collaboration, e.g., quotes a neutron excess of e = (N — Z)/ (N + Z) — 5.74% 
in their target, where N and Z are the (average) numbers of neutrons and protons in 
the target, respectively. The probabilities c p and c n to find a proton or neutron, are thus 
given by 

c P = \(l-e), c n = i(l + e). (5.4) 
The corresponding nuclear structure functions read 

fa(x) = c p a p (x) + c n a n (x), a = q,q, 7, (5.5) 

where a p (x) and a n (x) generically denote the PDFs for the proton and neutron, respec- 
tively. An isoscalar target obviously corresponds to the choice e = 0, i.e. c p = c n = \. In 
the following we make use of nuclear structure functions derived with the value for the 
neutron excess e quoted by NuTeV. 

6 Discussion of results 



6.1 Input parameters and setup of the calculation 

For the numerical evaluation we use the following set of SM parameters, 





= 1.16637 x l(T 5 GeV- 2 , 


a(0) 


= 1/137.03599911, 


a(M z ) 


= 1/128.93, 


M w 


= 80.425 GeV, 


M z 


= 91.1876 GeV, 


a s 


= 0.1172, 


M H 


= 115 GeV, 










m c 


= 0.51099892 Me V, 


«v 


= 105.658369 Me V, 


m T 


= 1.77699 GeV, 


m u 


= 66MeV, 


m c 


= 1.2 GeV, 


m t 


= 178 GeV, 


m d 


= 66MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV, 



(6.1) 

which essentially follows Ref. [30]. For the top-quark mass m t we have taken the value 
of Ref. [31]. The masses of the light quarks are adjusted to reproduce the hadronic 
contribution to the photonic vacuum polarization of Ref. [32]. 

As explained in detail in Section 5.1, we consistently use the set of MRST2004QED 
PDFs introduced in Ref. [10], which include 0(a) corrections. In order to quantify 
the modifications induced by these corrections to the PDFs we use the alternative set 
MRST2004QEDx, which we got from the MRST collaboration and which differs from 
the original MRST2004QED set only by neglecting the 0(a) effects. Since we do not 
include NLO QCD corrections in our calculation, it does not make sense to discuss the 
factorization scale dependence of our results, because the QED and QCD factorization 
scales are set equal in the MRST2004QED PDFs. As already done in Ref. [15], we set 
the QED factorization scale M 2 to the momentum transfer of the leptonic line; the only 
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difference comes from the fact that the MRST2004QED PDFs are not delivered for M 
values below their start scale M . Therefore, we use 

M 2 = max {-( Pl - hf, M 2 } , M 2 = 1.25 GeV 2 , (6.2) 

where pi and k\ are the incoming neutrino momentum and the outgoing neutrino or lepton 
momentum, respectively. Moreover, the PDFs require x > 10 -5 ; otherwise they are set 
to zero. 

The energy of the incoming neutrinos is set to E^ AB = 80 GeV, but the energy of all 
outgoing hadrons and photons is constrained by ^had+phot > 20 GeV, which differs from 
the setup used in Ref. [15]. 

As done in Refs. [14,15], we quantify the corrections to the R v ratio by 

^nc = — , oEcc = — — — , (6.3) 



where Scr^c/cc are the corrections to the NC and CC cross sections. These corrections, to 
a reasonable approximation, translate into a shift in the extracted on-shell weak mixing 
angle sin 2 6>w according to [14,15] 

AQin 2 fl - I - Sin2 9w + j Sin4 9w (SP V -J-AR" ^ (RA\ 
A sin U w - 40 , 2 (Ort NC + 0R CC ) . (6.4) 

27 w 

Technically we have extended our 0(a) calculation w.r.t. its description in Ref. [15] 
by implementing a generalization of the dipole subtraction formalism to non-collinear-safe 
observables. 1 All numerical results below are obtained with this technique, which turns 
out to be more efficient, but they have been confirmed by phase-space slicing as described 
in Ref. [15]. 

6.2 Numerical results on integrated cross sections 

Table 1 summarizes our numerical results as obtained in the three different IPS "a(0)" , 
"a(M z )" , and "G^" . In the first row, denoted u O(a), PDFqed" , we give the cross-section 
ratio Rq in lowest order as well as the corrections 5R^ C , SRq C , Asin 2 #w based on 0(a) 
corrections only. The results are based on the (9(a)-corrected PDFs MRST2004QED [10], 
but contributions from the 7-induced processes are not taken into account in this row. 
This means, apart from the different kinematical setup, the new PDFs (with DIS-like QED 
factorization), and the new input parameters, the first row is based on the calculation 
of 0(a) corrections described in Ref. [15]. The second row, denoted u O(a), PDFqedx", 
contains the corresponding results obtained with the set of PDFs MRST2004QEDx that 
differ from MRST2004QED only in the consistent neglect of the 0(a) corrections in the 
PDFs. Comparing the first two rows for each IPS, the influence of the implicit 0(a) 
corrections in the PDFs can be quantified to about 0.0004 in sin 2 W . All the remaining 
rows in Table 1 are calculated with the C(a)-corrected PDFs MRST2004QED. 

lr The dipole subtraction formalism as described for real QCD corrections [40] and photon radiation 
[34,35] is originally set up for collinear-safe observables where all final-state mass singularities cancel. 
The generalization to non-collinear-safe observables is partially described (for charged final-state particles 
only) in Ref. [36] ; the complete generalization will be published elsewhere [37] . 



11 



IPS 


contribution 


Rn 
U 


IN Lj 


5Rrr 


A sin 2 9 W 


a(0) 


£>(a),PDFqed 


0.30455 


0.0562 


-0.0840 


-0.0130 




0(a), PDFqedx 


0.31167 


0.0555 


-0.0841 


-0.0134 




univ. 2-loop 




0.01275 


-0.01279 


-0.00002 




univ. 2-loop' 




0.01229 


-0.01235 


-0.00003 




h.o. FSR 







0.000014(1) 


0.000007(1) 




h.o. FSR' 







0.000010 


0.000005 




7 process 




-0.000005 


-0.000063 


-0.000032 




"best" 




0.0690 


-0.0968 


-0.0130 


a(M z ) 


£>(a),PDFqed 


0.30455 


-0.0660 


0.0365 


-0.0138 




0(a), PDFqedx 


0.31167 


-0.0668 


0.0365 


-0.0142 




univ. 2-loop 




0.01248 


-0.01319 


-0.00033 




univ. 2-loop' 




0.01258 


-0.01330 


-0.00033 




h.o. FSR 







0.00073(1) 


0.00034(1) 




h.o. FSR' 







0.00072 


0.00034 




7 process 




-0.000005 


-0.000063 


-0.000032 




"best" 




-0.0535 


0.0240 


-0.0138 


Ga 


£>(a),PDFqed 


0.30455 


0.0016 


-0.0303 


-0.0135 




(9(a), PDFqedx 


0.31167 


0.0009 


-0.0303 


-0.0138 




univ. 2-loop 




-0.00039 





-0.00018 




univ. 2-loop' 




-0.00039 





-0.00018 




h.o. FSR 







0.00040(1) 


0.00019(1) 




h.o. FSR' 







0.00039 


0.00018 




7 process 




-0.000005 


-0.000063 


-0.000032 




"best" 




0.0012 


-0.0299 


-0.0135 



Table 1: Results on the ratio R v with ££ AB = 80 GeV, ^ a d+ P hot > 20 GeV in lowest order 
(Rq), corrections from NC and CC cross sections (5R^ C and 8Rq ), and shift Asin 2 6* w , 
for various input-parameter schemes (IPS). All but the "PDFqedx" numbers are obtained 
with the O (a) -corrected PDFs MRST2004QED; in PDFqedx the 0(a) corrections within 
the PDFs are switched off. The different rows correspond to the pure u O(a)" corrections, 
to improvements by "universal 2-loop" corrections, to "higher-order FSR" corrections, 
and to contributions from the "7 process" with a photon in the initial state. The number 
in parentheses indicates the statistical integration error in the last digit. For more details, 
see main text. 
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IPS 


contribution £ 


8R» CC 


A sin #vv 


a(0) 


h.o. FSR' 


0.3 


-0.000008 


-0.000004 






1 


0.000010 


0.000005 






3 


0.000040 


0.000019 


a(M z ) 


h.o. FSR' 


0.3 


0.00046 


0.00021 






1 


0.00072 


0.00034 






3 


0.00099 


0.00046 




h.o. FSR' 


0.3 


0.00024 


0.00011 






1 


0.00039 


0.00018 






3 


0.00055 


0.00026 



Table 2: Scale dependence of higher-order FSR corrections. 

The rows denoted "univ. 2-loop" and "univ. 2-loop' " show the impact of the universal 
electroweak corrections beyond 0(a) induced by Aa and Ap, as worked out in Section 3. 
In detail, the first of these rows correspond to the resummed versions of the two-loop 
effects [first equations in (3.9)— (3.11) and (3.14)-(3.16)], while the second corresponds to 
the genuine two-loop parts [second equations in (3.9) etc.]. The difference between the two 
versions is < 10~ 5 in sin 2 # w and, thus, negligible. Note that the size of these universal 
two- loop corrections to the cross sections is rather large, more than 1%, in the a(0) and 
a(Mz) schemes. In the G M scheme the universal effects beyond 0(a) are very small also 
for the individual cross sections, reflecting the fact that they are widely absorbed in the 
lowest order. Since the effects enter CC and NC cross sections in a similar way, they 
widely cancel in the cross-section ratio and influence sin 2 # w only at the level of 0.0003, 
depending on the IPS. 

The result on higher-order FSR corrections to the CC cross section are indicated by 
"h.o. FSR" and "h.o. FSR"' in Table 1, where the first of these rows result from the 
structure function as given in (4.2) and the second from an expansion of (4.2) up to 
0(a 3 ). The differences between these numbers, which quantify leading logarithmic FSR 
effects beyond 0(a 3 ), are negligible. The corrections in the a(Mz)- and G^-schemes are 
significantly larger than in the a(0)-scheme, because in those schemes the higher-order 
FSR corrections also receive a contribution from the 0(a) part, more precisely from the 
0(a) part calculated with a(0) minus the 0{a) part calculated with the a from the chosen 
IPS. In fact almost the entire effect is due to this change from a to ct(0) in the 0(a) FSR, 
which is 1.14% in SRq C when calculated with a(0). In the scheme, higher-order FSR 
affects sin 2 9 W at the level of 0.0002. 

As explained in Section 4, our results depend on the QED splitting scale Q for collinear 
FSR in 0(a 2 ) and beyond, while in 0(a) the scale Q is fixed by the non- logarithmic 
corrections. Table 2 illustrates the residual scale dependence by varying the parameter 
£ in (4.6) between 0.3 and 3; the default choice is £ = 1. The table shows that in the 
a(0)-scheme the scale uncertainty is of the order of 0.00002 when translated into a shift in 
sin 2 6*w- While this number corresponds to genuine FSR effects beyond 0(a), the larger 
scale dependence in the a(M z )- and G M -schemes originates from the scale change in the 
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source 


A sin 2 #w 


0(a) effect in PDFs 


0.0004 


univ. electroweak effects 


0.00003 


h.o. FSR 


0.00002 


non-univ. electroweak effects 


0.0003 



Table 3: Estimates of theoretical uncertainties from missing electroweak higher-order 
effects in the G M -scheme. 

term [a(0)/n — a/ir] ln(m 2 /Q 2 ), which changes the coupling a to a(0) in the 0(a) FSR 
corrections. 

The contributions of the 7-induced processes are shown in Table 1 in the rows "7 pro- 
cess". Since we have set the photon coupling to a(0) in all three IPS, the relative impact 
of these corrections is identical in all the IPS. The effect on sin 2 # w is about 0.00003 and, 
thus, negligible. 

The numbers called "best" in Table 1 show the sum of the O(a) and higher-order 
corrections, where the latter are the sum of the contributions labelled "univ. 2-loop", 
"h.o. FSR"', and "7 process". The fact that these "best" numbers agree with the purely 
C(o;)-corrected results up to < 0.0001 in sin 2 # w is accidental, since universal effects and 
FSR, which widely cancel each other, are of different physical origin. This means, in 
particular, the cancellation will not take place for a different kinematical setup (neutrino 
energy, hadronic energy cut, etc.). 

Finally, in Table 3 we estimate the theoretical uncertainties from missing electroweak 
higher-order effects from the results discussed above. Table 1 clearly shows that the 0(a) 
corrections implicit in the PDFs affects sin 2 # w of the order 0.0004, which is not negligible 
w.r.t. the error of about 0.0016 quoted by NuTeV. (9(a)-corrected PDFs should therefore 
be used in such an analysis. If this is done, the error estimate of 0.0004 will of course 
be reduced. In Table 3 we estimate the missing universal electroweak effects in the G^- 
scheme, which is the most reliable IPS, by changing the parameter Ap by 2% of its size, 
which is of the order of effects beyond two loops [38]. The resulting uncertainty is about 
0.00003 and thus negligible. Missing higher-order effects from FSR have been quantified 
to be of the order 0.00002 in Table 2 above. We take this value from the a(0)-scheme, 
because the larger uncertainties from changing a to a(0) in the two other schemes are part 
of non-universal electroweak effects beyond 0(a). These non-universal electroweak effects 
are estimated as follows. In the G M -scheme, the 0(a) corrections to the R v ratio are about 
2.9% where about 1.1% are due to FSR (£ = 1). The relative coupling of FSR is clearly 
ct(0), but the effective value of a in the non-universal 1.8% can only be fixed in a two- 
loop calculation. Therefore, changing the value of a^^ to a(0) or a(Mz), i.e. varying the 
1.8% by about 3%, should give an idea about missing electroweak two- loop effects. This 
procedure changes sin 2 6*w by 0.0003, which means that missing non-universal electroweak 
effects should dominate the theoretical uncertainty from electroweak corrections. 
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6.3 Numerical results on differential cross sections 

Figures 1-3 illustrate the differential cross sections da/dy and d 2 a/dx/dy and the 
corresponding corrections. The curves labelled "total" correspond to da/dy shown as 
histograms, the other curves correspond to d 2 a/dx/dy integrated over the indicated x 
range. Note that we actually distinguish between the variable x, which is the argument 
of the PDFs, and the Bjorken variable x-Bj, which is reconstructed from initial- and final- 
state variables. The lowest-order predictions of Figure 1 show that the integrated cross 
sections, both for the NC and CC channel, are dominated by the x range (x = xbj in 
lowest order) x < 0.4, while the distributions in y are flat. The corrections to the NC 
process (Figure 2) vary between —0.1% and +0.5% for all x, XBj, and y values. The 
distributions in x and x-q] are almost identical, and the corrections in the various x and 

ranges do not deviate by more than 0.2% from the correction to the corresponding 
fully integrated quantity da/dy. The situation is very different in the CC shown 
in Figure 3. The corrections to da/dy vary within ~ 2—4%. The upper left plot in 
Figure 3 demonstrates that the corrections vary very little in x, but the upper right plot 
shows a very large variation in xb v This is due to the fact that binning w.r.t. the parton 
x does not destroy the inclusiveness of collinear FSR off the muon, but binning in x B j 
does. In other words, a cut in x does not cut into the cone of collinear FSR, but a cut in 
rr B j does. These features also show up in higher-order FSR as depicted in the lower plots 
of Figure 3. 

The results on the distributions are particularly interesting for the issue of implement- 
ing electroweak corrections in Monte Carlo generators upon reweigthing certain distribu- 
tions. The corrections to the differential CC cross sections suggest that a reweighting in 
distributions w.r.t. x is less delicate than in x-Qy 

7 Conclusions 

A thorough understanding of electroweak corrections to deep-inelastic neutrino scat- 
tering is indispensable for exploiting measurements of NuTeV and NOMAD to extract 
a value for the on-shell weak mixing angle sin 2 #w = 1 — M^/M^- In this paper, we 
have supplemented our previous calculation of electroweak 0(a) corrections by higher- 
order effects, in order to reduce the theoretical uncertainty from unknown electroweak 
corrections. 

We take into account universal two- loop effects from Act and Ap as well as higher-order 
final-state photon radiation off muons in the structure function approach. The impact 
on sin 2 # w of these effects are about 0.0003; compared to the experimental uncertainty 
quoted by NuTeV, which is about 0.0016, these effects are small. 

Moreover, we employ the recently released (9(o;)-improved parton distributions 
MRST2004QED and identify the relevant QED factorization scheme, which is DIS like. 
Switching on and off the implicit 0(a) corrections in the parton distributions, changes 
sin 2 6 W at the level of 0.0004. 

As a technical byproduct, we describe two variants of phase-space slicing and the dipole 
subtraction method for calculating processes with initial-state photons. Such photon- 
induced processes become relevant, because the 0(a) evolution of the parton distributions 



15 



dy 



[fb] 



10 - 



1 - 



0.1 r 



total 



0.KK0.2 
0.0<x<0.1 
0.2<x<0.3 

■0.3<x<0.4 

0.4<x<0.5 

0.5<z<0.6 
0.6<z<0.7 



_L 



0.2 



0.4 0.6 

y 



NC 



1 1 

total 


1 1 

CC " 


- 0.KK0.2 
















1 1 


1 1 



1 



0.2 



0.4 0.6 

y 



0.8 



Figure 1: Lowest-order predictions for the distributions in y (total) and individual contri- 
butions from different ranges in x (= XBj in lowest order). The NC channel for v scattering 
is shown on the l.h.s., the CC channel on the r.h.s. 
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Figure 2: Relative electroweak corrections 5 U = 5a u /<Tq to the distributions in y (total) 
and to individual contributions from different ranges in x (l.h.s.) and in XBj (r.h.s.). 
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Figure 3: Relative electroweak corrections 5 U = 5a v /ctq (upper plots) and higher-order 
FSR (lower plots) to the distributions in y (total) and to individual contributions from 
different ranges in x (l.h.s.) and in x-q s (r.h.s.). 
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generates a photon distribution. Owing to the smallness of the available scattering energy, 
the impact of photon-induced processes on sin 2 9^ is negligible for NuTeV. 

A numerical discussion of the higher-order effects suggests that the remaining theo- 
retical uncertainty from unknown electroweak corrections is dominated by unknown non- 
universal two-loop effects and is of the order 0.0003 when translated into a shift in sin 2 W . 
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Appendix 

Calculation of photon-induced cross sections 

As explained in Section 5.1 in detail, a consistent evaluation of 0(a) corrections to 
quark-induced hadronic processes includes an evaluation of related photon-induced chan- 
nels. In neutrino DIS we have the leading-order processes v^q — > v^q and u^q — > fi~q', 
with q and q' generically denoting the relevant quarks and antiquarks. The corresponding 
photon-induced channels are — > v^qq and — > pTq'q. The matrix elements for these 
reactions obviously follow upon crossing from the bremsstrahlung processes v^q — > v^q^ 
and v^q — > /J,~q ,r y, which already entered the calculation of 0(a) corrections described 
in Ref. [15]. Alternatively they can also be derived upon crossing the results for the real 
0(a) corrections to Drell-Yan-like W production as explicitly given in Ref. [39]. We 
set the coupling factors such that the squared matrix element is proportional to a 2 o;(0), 
where a is equal to cc(0), at(Mz), or ac^ according to the respective IPS; the factor a(0) is 
taken independent of the IPS, because this is the correct effective coupling for the incom- 
ing photon which is on shell. Since the matrix elements have been derived for massless 
fermions the integration over the domains of collinear initial-state splittings 7 — > ff re- 
quires a restoration of finite-fermion-mass effects. The resulting mass-singular logarithms 
of initial-state quarks are absorbed into the photon PDF as described in Section 5.1. 

In the following we describe three simple methods to restore such effects in a generic 
reaction a{p a ) +j(k) — > f(pf) + A, where a stands for any massless incoming parton and / 
is the light fermion or antifermion whose mass effects we are interested in. The momenta 
of the particles are given in the parentheses. The remainder X may contain additional 
light fermions which can be treated in the same way as /. The results obtained with the 
three methods are in mutual numerical agreement, where the dipole subtraction method 
yields the smallest integration error when using the same statistics in all the methods. 

(i) Effective collinear factor 

The collinear singularity in the squared matrix element \A4 a7 ^fx\ 2 occurs if the angle 
9f between / and the incoming 7 becomes small; in this limit the scalar product (kpf) 
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is of 0(mj), where rrif is the small mass of /. Neglecting terms that vanish in the limit 
rrif — > the polarization-summed squared matrix element (including the average over 
initial-state polarizations) asymptotically behaves like [40] 



£1^ 



QW 



aj^fX 



(Pa,k)\' 



r 



pol 



k Pf^° x(kpf) 



1 - 2x(l -x) + 



m 



kp f 



x 



J2\M af -^ X (p a ,Pf = xk)\ 2 , 
pol 

(A.l) 

where x — l—p®/k°, Qfeis the electric charge of /, and M. a j^x * s the matrix element of 
the related process which results from a-y{— > aff) — > fX upon cutting the / line if this 
process proceeds in the collinear limit mainly via the 7 — > // splitting. The incoming 
momenta relevant in the different matrix elements are given in parentheses. 

The factorization formula (A.l) can be used to relate the two versions for the squared 
matrix element \M. a ~f-*fx\ 
eliminating \M. a f^ x \ 



£ \M ai ^ f x 

pol 



2 obtained with a finite or vanishing fermion mass rrif, just by 

(A.2) 



1 The result is 

(Pa, k)\ 2 



9(m f ,x,k°,e f ) J2 \M a ^ f x(Pa,k)\ : 

pol 



mf=Q 



with the "effective collinear factor" 
g(m f ,x, k°, Of) 



{l-xf 4(A;°) 2 sin 2 ^ 


m 2 + 4(A;°) 2 sin 2 ^(l- 


- x) 2 ((1 - x) 2 + x 2 ) 


(1 — x) 2 + X 2 


m 2 +4(A;°) 2 sin 2 ^(l 


-x) 2 


2 



(A.3) 



The squared matrix element on the l.h.s. of Eq. (A.2) can be integrated over the phase 
space for a massless /, yielding the correct result for a massive / up to terms that are 
suppressed by powers of rrif. This means that the collinear singularity is regularized by 
the correct cutoff term In rrif supplemented with the correct constant contributions for 
rrif — > 0. However, there are spurios terms of 0{rrif/Q) where Q is a typical scale in the 
process. 



(ii) Phase-space slicing 

In the collinear limit 9f 
can be factorized as 



/ d$ 

Je f <A9 



aj^fX 



the phase-space volume element d$ a7 _>/x(Po, k) 

(l-x)(k ) 2 rA6 



0, rrif 

,, I dx J**af^x(Pa,xk) y2 



(P°> k ) A^O L dX I d ®af^x(Pa,xk) K± jf d^SHl^, 

(A.4) 

where A6 <C 1 is a small cutoff parameter defining the collinear region. Here we also used 
that the incoming parton a is massless. Making use of this factorized phase space and of 
the asymptotic form (A.l) of the squared matrix element, the integral over small emission 
angles 9f can be carried out, yielding 



fx{Pa, k) ^ <Ae ^ £ 

So -<Ae d *«^/*fo»' ^ £ \ M ^fx(Pa, k)\ 
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QfCH f 1 n 1 f 

J o dx H(m f , x, k , AO) J d$ af ^ x (p a , xk) 

xY.\Maf^x(Pa,Pf=xk)\ 2 + O(A0), 



2ir 



pol 
)2, 



NcJ ^B j ldxH ( m f^^^ Ad )^x(Pa,Pf = xk) + O(A0) 

(A.5) 



with the auxiliary function 

i2/l0\2/1 _ „.\2 



/ Af) 2 (k°) 2 ('\ — r) 2 \ 
H(m f ,x, k°,A0)=P fl {x) \n[ '-)+ 2x{l - x) 



(A.6) 



containing the collinear singularity. The splitting function Pf 7 (x) is defined in Eq. (5.2), 
and N c j is the colour factor for the fermion / (A^ C) i epton = 1, iV C)quark = 3). 

The remaining phase-space integration over the non-collinear region Of > AO can be 
carried out for a massless / (provided no other singularities are connected with /) without 
modifications. 

(Hi) Dipole subtraction 

The general idea of a subtraction method is to subtract and to add a simple auxiliary 
function from the singular integrand. This auxiliary function has to be chosen such that 
it cancels all singularities of the original integrand so that the phase-space integration of 
the difference can be performed numerically. Moreover, the auxiliary function has to be 
simple enough so that it can be integrated over the singular regions analytically, when 
the subtracted contribution is added again. The dipole subtraction method, originally 
worked out for massless QCD [33], provides a general algorithm for the construction of 
this auxiliary function and for its integrated counterpart. 

The dipole subtraction method for photon radiation off massive and massless fermions 
has been described in Ref. [34] and Ref. [35]. Following the same strategy, we work out a 
similar formalism for the collinear splitting 7 — > // of a photon in the initial state into a 
light // pair. More generalizations of the dipole formalism of Ref. [34] will be published 
elsewhere [37]. 

The function that is subtracted from the integrand ^ po i \M a y^fx(pa, k)\ 2 , which is 
derived for rrif = 0, is defined as follows, 2 

|^ sub | 2 = Q)e 2 h ja (k,p f ,p a ) \M af -^ x (Pa,Pf = x ia k, k x )\ 2 , (A.7) 

pol 

with the radiator function 

7 /, x P fj( X ia) Pak-pfk-p a pf 

hja(k,Pf,Pa) = , x ia = J — J -. A.8 

XjaiPfk) p a k 



2 This construction closely follows the case of an "initial-state emitter" and of an "initial-state specta- 
tor" in the case of photon or gluon radiation, as described in Refs. [33-35]. 
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The final-state momenta entering the squared matrix element on the r.h.s. of Eq. (A. 7), 
generically denoted kx, follow from the original momenta k x upon the Lorentz transfor- 
mation 

k x = A"„A& (A.9) 

with 

( P 4- PW P 4- P) 2P^P 

(A.10) 

It is straightforward to check that |-M S ub| 2 possesses the same asymptotic behaviour as 
Epoi \M ai ^fx\ 2 in Eq. (A.l) for m f = 0. Thus, the difference £ pol \M ai ->fx\ 2 - |-M sub | 2 
can be integrated numerically for nif = 0. The correct dependence of |.M su b| 2 (and the 
related kinematics) on a finite mj is, however, needed when this function is integrated 
over Of leading to the collinear singularity for Of — > 0. The details of this part of the 
calculation will be described elsewhere [37]; here we give only the final result: 

b Q 2 f a f 1 

a™^ fx (p a ,k) = N cJ ^— &xH{mf,x,p a k)a a f^x{Pa,xk) (A.ll) 

with the auxiliary function 

H(m f ,x,p a k) = P h {x) \n( 2Pak ^~ x)2 )+2x{l-x). (A.12) 

V 171 f J 



The contribution a s ^t^f X simply has to be added to the result for the cross section obtained 
from the integral of the difference J2 po \ \M.a^fx\ 2 — |-M S ub| 2 - 
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